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Abstract — Using the EXIT chart approach, a necessary and 
sufficient condition is developed for the local stability of iterative 
decoding of multi-edge type (MET) doubly-generalized low- 
density parity-check (D-GLDPC) code ensembles. In such code 
ensembles, the use of arbitrary linear block codes as component 
codes is combined with the further design of local Tanner graph 
connectivity through the use of multiple edge types. The stability 
condition for these code ensembles is shown to be succinctly 
described in terms of the value of the spectral radius of an 
appropriately defined polynomial matrix. 



I. Introduction 

Multi-edge type (MET) low-density parity-check (LDPC) 
codes were originally proposed in [1] as a framework to 
capture both degree- 1 variable nodes (VNs) and punctured 
bits in the analysis of LDPC code ensembles, and to achieve a 
finer control of the connectivity between VNs and check nodes 
(CNs) in the ensemble definition. The new framework allowed 
the design of powerful finite-length LDPC codes over the 
additive white Gaussian noise (AWGN) channel, with a very 
good compromise between waterfall performance and error- 
floor. Since then, several aspects of MET LDPC codes have 
been investigated such as, for instance, their average weight 
distribution [2]. Traditional unstructured irregular LDPC code 
ensembles parametrized through their degree distribution pair 
[3] may be seen as MET ensembles where all edges in the 
Tanner graph are of the same type. On the other hand, LDPC 
codes based on protographs [4] may be seen as MET LDPC 
codes such that no two edges connected to the same VN or 
to the same CN are of the same type. 

Another way to extend the original framework of unstruc- 
tured LDPC code ensembles consists of replacing the VNs and 
the CNs with linear block codes other than repetition codes and 
single parity-check (SPC) codes, respectively. The resulting 
LDPC-like codes are called doubly-generalized LDPC (D- 
GLDPC) codes [5], and extend the original idea of generalized 
LDPC (GLDPC) codes [6], where only the CNs were replaced 
with generic linear block codes. Several theoretical aspects 
of unstructured D-GLDPC codes have been recently clarified, 
such as their stability condition over the binary erasure channel 
(BEC) [7], and the analysis of the exponent of their weight 
distribution [8]. 
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The two different extensions can be considered together, 
leading to the concept of MET D-GLDPC code ensemble. This 
represents a very general framework for design and analysis 
of LDPC-like codes, enabling to handle different variable and 
check component codes along with puncturing; VNs and CNs 
with local minimum distance 1, including degree-1 VNs (state 
variables), can be also considered. The asymptotic weight 
enumerators for MET D-GLDPC codes were investigated in 
[9], while EXIT analysis to calculate the threshold of MET 
D-GLDPC codes over the BEC was developed in [10]. 

In this paper, we analyze the convergence properties of 
the belief -propagation (BP) decoder for MET D-GLDPC code 
ensembles over the BEC by developing its stability condition, 
i.e., the condition under which the erasure-free state attracts 
the decoder, in the asymptotic setting where the codeword 
length tends to infinity. If and only if the condition is sat- 
isfied, BP decoding can in principle succeed, provided the 
BEC erasure probability is below the threshold which can be 
calculated using the technique in [10]. The stability condition 
is obtained in the case where the are no punctured bits and 
where the local minimum distance of each VN and CN is at 
least 2. It can be shown that the obtained condition coincides 
with that developed in [1], [11, Ch. 7] in the special case of 
MET LDPC codes. 

II. Preliminary Definitions 

A. Concept of D-GLDPC Codes 

A D-GLDPC code consists of a set of CNs and a set of VNs. 
Each CN corresponds to some arbitrary linear 'local' code. On 
the other hand, each VN corresponds to some arbitrary linear 
'local' code, together with its encoder (i.e., generator matrix). 
Graphically, each CN and each VN may be viewed as having a 
set of sockets corresponding to the bits in the local codeword. 
The sockets of the CNs are connected by edges to the sockets 
of the VNs in a one-to-one fashion; the resulting graph is 
called the Tanner graph of the D-GLDPC code. A codeword of 
the D-GLDPC code is an assignment of an information word 
to each VN such that the local encoding of this word at each 
VN assigns an encoded bit to each edge of the Tanner graph 
in such a way that the resulting configuration forms a valid 
local codeword from the perspective of every CN. It is easily 
seen that if the local code at each CN is a single parity-check 
code and if the local code at each VN is a repetition code, the 
resulting D-GLDPC code reduces to an ordinary LDPC code. 

B. MET D-GLDPC Code Ensemble Definition 

In MET D-GLDPC codes, we distinguish between rig 
different edge types. Each edge type is identified by an index 



in the set S' ~ {1, 2, . . . , ng}. Furthermore, we distinguish 
between different VN types and different CN types. Each VN 
type is identified by a triplet 7 ~ (u-y, 6^, d^), where: 

> identifies a {q^,k^) variable component code, where 
is the code length and the code dimension, and a 
specific encoder (i.e., generator matrix G-^) of it; 

• b-y is a binary vector of length which specifies the 
local puncturing pattern for the VN. Specifically, for 
i S {1, 2, . . . , fc^}, if b^,i = then the corresponding 
encoded bit of the D-GLDPC code is punctured, and it 
is not punctured otherwise; 

• is a vector of length whose i-th element d^,i G S" 
specifies the edge type of the i-th VN socket. 

Each CN type is identified by a pair 5 = {cs,ds), where: 

• cs identifies an {ss,hs) check component code (regard- 
less of its representation), where ss is the code length 
and hs the code dimension; 

• is a vector of length ss whose i-th element ds^i G S" 
specifies the edge type of the i-th CN socket. 

The set of all VN types 7 is denoted by ,!^v, and set of 
all CN types S by ,!^c- Moreover, the fraction of edges of 
type I £ connected to VNs of type 7 £ is denoted 
by A^,;, while the fraction of edges of type I G S connected 
to CNs of type 5 G ,9^0 by ps.i- We have A^,; > if and 
only if the generic VN of type 7 has at least one socket of 
type I, and A^,; ~ otherwise. An analogous statement can 
be made regarding p^j. Also, q^,i and ss,i denote the number 
of sockets of type / for a VN of type 7 and for a CN of type 
(5, respectively. The constraints 'Yliies ~ 97 ^7 ^ '^v and 

Y^ies ^^.i = •'^'5 ^ ■^c hold. 

An example of D-GLDPC code ensemble is depicted in 
Fig. [T] As opposed to single-edge type codes, where a unique 
edge interleaver for all edges is present, for MET codes a 
dedicated edge interleaver is present for all edges of the 
same type. For a given codeword length, each code in the 
ensemble corresponds to a specific realization of the Tig edge 
interleavers, where all realizations of each edge interleaver are 
equiprobable. 

In the following, we make the assumption that there are no 
VNs and CNs with local minimum distance 1. We also make 
the assumption that no encoded bit of the D-GLDPC code 
is punctured, i.e., for each VN type 7 = {v^,b^,d^) G J^v 
the vector has no '0' entries. Finally, we assume that no 
variable or check component code has idle bits. 

C. Further Definitions 

Throughout the paper, vectors are intended as column 
vectors. We define and 1 as the vectors of length rip whose 
elements are all equal to and all equal to 1, respectively. 
Morever, we define le as the vector of length rig whose 
elements are all equal to except the element in position 
e which is equal to 1. The subset of VN types 7 with local 
minimum distance 2 is denoted by C and the subset 
of CN types 5 with local minimum distance 2 by Q ^c- 

For b G =^C2 and l,m G (a, we denote by £,2^\l,m) the 
number of ordered pairs of sockets of a CN of type S, such 
that the first socket is of type I and the second of type m, and 




Fig. 1 . Example of a MET D-GLDPC code ensemble. The set of edge types 
is ^ = {1, 2, 3}, where edges of type 1, 2 and 3 are depicted in red, blue and 
green, respectively. A separate edge-interleaver is present for each edge type. 
The set of VN types is = {7i,72i73}> where the VNs of type 71, 72 
and 73 are represented as cyan, brown and pink circles, respectively. The VNs 
of type 71 are (1, 1) codes, those of type 72 are (2, 1) codes, and those of 
type 73 are (3, 2) codes. Each VN of type 71 and 72 is associated with one 
encoded bit of the D-GLDPC code, each VN of type 73 with two encoded 
bits of the D-GLDPC code. The encoded bits corresponding to type-71 VNs 
are punctured. The set of CN types is = {<5i,i52}, where the CNs of 
type (5i and 62 are represented as yellow and magenta squares, respectively. 
If each CN of type 5i introduces two parity-check equations and each CN of 
type 82 one parity-check equation, then the overall code dimension is K = 8. 
Since four bits are punctured, the overall codeword length is A'^ = 28, and 
the code rate is ij = 2/7. 

such that the assignment of a '1' to these sockets and a '0' to 
all other CN sockets results in a (weight-2) local codeword. 
Note that ^2'^' (Z, m) = ^2'''' ("^j 0- for 1,iti £ r^, we define the 
nonnegative real parameter C''™ as 

If CNs of type S G ^C2 have no sockets of type / (in which 
case ps^i = and ssj = 0), then we set C*''™ = by 
definition. Denoting by A^\l,m) the number of weight-2 
local codewords of a CN of type 5 G =^C2 such that one of 
the two '1' local encoded bits corresponds to a socket of edge 
type I and the other to a socket of edge type to, we have 
^^2^\l,m) = 4*Hz,to) for / 7^ to and = 24''' (^,0- 

For 7 G ^V2 and l,m G S", we denote by X2'l!(^i"^) the 
number of ordered pairs of sockets of a VN of type 7 such 
that the first socket is of type / and the second of type m, and 
such that the assignment of a '1' to these sockets and a '0' 
to all other VN sockets results in a (weight-2) local codeword 
generated by a local input word of weight u. Similarly to the 
CN case, we have X2''2(^'™) ~ X2''2("^jO- The nonnegative 
polynomial P''™(a;) (with real coefficients) is defined as 

P^^-ix):^ (^)j:x^^,l{l,m)x". (2) 

If VNs of type 7 G ^V2 have no sockets of type I (in which 
case A.^^; = and q^^i = 0), then we set P''™(x) = by 
definition. Moreover, denoting by B':^l^{l,m) the number of 
weight-2 local codewords of a VN of type 7 G .'^V2 generated 



by local weight-w input words, and such that one of the two '1' 
local encoded bits corresponds to a socket of edge type I and 
the other to a socket of edge type m, we have X2'«(^i"^) = 
Bi-^lil m) for I ^ m and x^^}Xl, = 25^:2 (Z, 0- ' 

Remark 2.1 : For l,m G S" and Z 7^ to, in general we 
have C'-" ^ C""'' and P'^'"(a:) ^ P'"-'(a;). 

Remark 2.2: For single-edge type codes, <o' is a singleton 
fS" = {I}, P^'^{x) reduces to the polynomial P{x) and C''' to 
the parameter C characterizing ordinary D-GLDPC codes [7], 
[8]. For codes constructed from protographs, no two sockets 
of a VN or CN are of the same type. Hence, in this case we 
have P^'\x) = and C'-' = for all I G ^. 



D. Multi-Type Information Functions 

Although in this paper we make some assumptions on 
the VNs and CNs (see the last paragraph of Sec. III-Bb . the 
definitions provided in this subsection are more general and 
do not rely on such assumptions. 

Consider a CN of type 5 ~ {cs,ds) G and let Gg 
be any generator matrix for the associated component code. 
From G^, form matrices G^^;, where Ggj is the {hs x ss^i) 
matrix composed of the columns of Gs associated with the 
bit positions of type I G (irrespective of the order of these 
columns). Then, for any integer rig-tuple g = (gi , (72, . • . , gn^ ) 
satisfying < gi < ss,i for all I G (o, the g-th multi-type 
information function of the CN is defined as 



;(<5) 



;(<5) 
g 



(3) 



where Sg'^'' is a matrix formed by selecting gi columns in Gs,i 
(irrespective of the order of these columns) and where y^a(*) 

denotes the summation over all YVi=i (^g,'0 matrices Sg'^-'Q 

For a VN of type 7 = {v^, b-y, d^), let G-y be the specific 
generator matrix identified by v^. Moreover, let be the 
Hamming weight of b^. From G-y, form Tie matrices G-y,/, 
where Gj^i is the {k^ x q^ j) matrix composed of the columns 
of G^ associated with the bit positions of type I G S" 
(irrespective of the order of these columns). Then, for any 
integer rie-tuple g = (gi, 32, • ■ • , 5nJ satisfying < gi < q^.i 
for all I G (o, and for any integer < u < fc^, the (g; it)-th 
multi-type split information function of the VN is defined as 



g(7) .- 



rank 



ail) 

'-'g-.u 



(4) 



multi-type split information function of a type-7 VN depends 
on the local mapping between information and encoded bits 
defined by G^. It also depends on the local puncturing pattern 
defined by b^. 

III. Exit Analysis and BP Decoder Stability 

EXIT analysis of a MET D-GLDPC code ensemble with Ue 
edge types consists of modeling the average behavior of the 
iterative decoder, in the asymptotic case where the codeword 
length tends to infinity, through an ng-dimensional discrete 
dynamical system tracking Tig average extrinsic information 
values, one for each edge type. For e £ S", the e-th value we 
track is the average extrinsic information over the edges of 
type e, outgoing from the VN set towards the CN set. 

Let £ > I denote the decoding iteration index. Let /|;yg 
and g be the average extrinsic information over the edges 
of type e outgoing from the VN set and from the CN set, at 
the £-th decoding iteration, respectively. Moreover, let l\y ^ 
and g be the average a priori information over the edges 
of type e incoming towards the VN set and towards the CN 
set, at the £-th decoding iteration, respectively. EXIT analysis 
equations of a MET D-GLDPC code ensemble over a EEC 
with erasure probability e may be expressed as 



and 



^EC,1 



Ievi(/ 



'-1 



£-1 



— 'EVriel-'Ay,!' ^Ay,2' ■ ' ■ ' ^Ay.rie' 



(5) 



(6) 



^AC,l'-^AC,2' ■ ■ ■ ' ^AC,nJ 



The 2ne equations Q and (|6]l, together with /^^ ^ = I^y ^ 
Vi G S", I^^y^ = /Ic-.i Vi G S", and /^^^ = Vi G define 
a recursion that can be expressed in the compact form 



(7) 



for £ > and where Iev = [lEV.iJEV.2,---JEV,n,]'^ 
is a column vector whose elements are the rig values to be 
tracked. The jig -dimensional discrete dynamical system ^ 
models the asymptotic (in terms of codeword length) evolution 
of the BP decoder over a EEC with erasure probability e. 
The function /(•) can be evaluated exploiting results devel- 
oped in [10]. In more detail, neglecting the iteration index 
L for e G ^ we have Iev, e(-fAy,i j -fAy,2, ■•■ 7 -^Ayne , £) — 
E^g^v ^7,6 Iev^^^v-,!: ^Ay2, ■ • ■ , lAV,n, , e) where 



where SgTu is a matrix formed by selecting gi columns in G^./ 
(irrespective of the order of these columns) and u columns 
among the columns of I^^ (order-A:^ identity matrix) 
corresponding to the support of b^. In (01), X^s'"*' denotes 



the summation over all 



n: 



1=1 



7 ■ 

i-tA 



matrices S 



(7) 
g;M- 



While the g-th multi-type information function of a type-S 
CN is independent of the specific choice of Gs, the (g; M)-th 

'if for some I £ «f the CN has no sockets of type /, then gi is conventionally 
set to 0. This convention shall be adopted also for the multi-type split 
information functions defined for the VNs. 



(7) 
EV,, 



{Iav,i, Iav,2, ■ ■ ■ , Iavu^,^) = 1 



z=0 



X ^(1-/Ayi)*^(/Ayi)'-^-*^- 

ti=0 

X {l-lAV,ef'{lAV,eY-'--^- 



^ (1 - -^Ay.ne)*"" {lAV,n^) 



(8) 



and 

Moreover, we have IecA^ac,!, Iac,2, ■ ■ ■ , Iac^uJ = 
EiGJTc Ps,e \ecMaca,Iac,2, • • • , /^anj, where 

^EC,ei^AC,l, IaC,2, ■ ■ ■ , lAC,nJ 
ti=0 

X {l~lAC,eY'{lAC,eY'--^-'^--- 



te=0 



E (1 - ^ac,„J*-(/ac,„J^^ ^-*- ^) (10) 



and 



(11) 



Lemma 3.1: For a MET D-GLDPC code ensemble such 
that all VNs and CNs have local minimum distance at least 
2 and such that no encoded bit is punctured (6-y is the all-1 
vector for all 7 G ^v), [Iev,i,Iev,2^ ■ ■ ■,lEV,nJ'^ = 1 is a 
fixed point of O regardless of e, i.e., /(I, e) = 1 Ve G (0, 1)0 

The fixed point Iev = 1 corresponds to a state of 
the system in which no erasure messages are exchanged 
between the VN set and the CN set, i.e., in which all 
encoded bits of the D-GLDPC code are known. A trans- 
mission over the BEC may be then modeled as a pertur- 
bation of the system state from Iev = 1 to Iev = 
[lEV,i(0,0,...,0,e),...,lEV,„,(0,0,...,0,e)]^ I%y{e), 
and the corresponding evolution of the BP decoder is modeled 
by (|7]i. Decoding is successful when, starting from J^y(e), we 



have liuif 



'By 



1. In order for the limit to be 1, it is 



necessary that the steady-state equilibrium Iev = 1 acts as a 
local attractor for the system, or, equivalently, that it is locally 
stable. The stability condition is established in the following 
theorem, which represents the main contribution of this paper. 

Theorem 3.1: Consider a MET D-GLDPC code ensem- 
ble with rie edge types. Assume that there are no VNs and 
CNs with local minimum distance 1 and that no encoded bit 
of the D-GLDPC code is punctured (6^ is the all-1 vector for 
all 7 G ^v)- Define C as the (rie x rig) nonnegative matrix 
whose {I, m)-th entry is (7''™ in ([T]i. Moreover, define P{e) as 
the (rie X rie) nonnegative matrix of polynomials whose [l, m)- 
th entry is P''™(e) in Q. Then, the fixed point I ev = 1 of 
O is locally stable if and only if 



a(P(e)C)<l 



(12) 



^The proof of Lemma [xT] is omitted due to space constraints. Tlie lemma 
can be easily proved by proving that, in Js) and UQK under the two mentioned 
hypotheses, 1^^^^^ -i- 1 as {lAV,i,lAV,2,---jAV,n^) (1,1, ■•■,!) 
V(7,e) e ,J?v X and l^^^'e ^ ^ ''^ (/ac,i, ^AC,2, • • • , -f^CnJ ^ 
(1,1, V(<5,e) G .^c x>. 



being cr{A) the spectral radius of a square matrix A, i.e., the 
magnitude of the eigenvalue of A with the largest magnitude. 

Interestingly, inequality (fT2] i represents the "natural" exten- 
sion to the MET framework of the condition P{e)C < 1 
proved in [7] for the single-edge type case. A sketch of proof 
of Theorem [TT] is provided in Section HV] Theorem l3 . 1 1 allows 
us to develop a simple sufficient condition for local stabihty 
of fixed point Iev = 1, as follows. 

Corollary 3.1: Consider a MET D-GLDPC code ensem- 
ble with rip edge types. Assume that there are no VNs and 
CNs with local minimum distance 1 and that no D-GLDPC 
encoded bit is punctured. Moreover, assume that the follwing 
condition is satisfied: If a socket of VN of type 7 G ^V2, 
associated with the support of a weight-2 local codeword, is 
of type I G (f, then for all S G -^C2 a CN of type S has 
no sockets of type I associated with the support of a weight- 
2 local codeword. Then, the fixed point Iev = 1 of dTji is 
locally stable for any BEC erasure probability e. 

Proof: Simply observe that in this case P{e) C is the 
all-zero matrix and, consequently, cr{P{€) C) = for all e. □ 

IV. Sketch of Proof of Theorem I3. II 

For ease of presentation, we consider the case = 2. The 
extension of the proof to the case of ?ie > 2 edge-types is 
straightforward. 

It is well-known that the local stability of a fixed point 
of a multidimensional discrete dynamical system such as ^ 
depends on the eigenvalues of the Jacobian matrix calculated 
in the fixed point. Specifically, the fixed point is a local 
attractor when the magnitude of all eigenvalues of the Jacobian 
matrix is less than 1 or, equivalently, if and only if the spectral 
radius of the Jacobian matrix is less than 1. Hence, we need to 
prove that J(l, e) = -P(e) C, where J(l, e) is the {ue x rie) 
Jacobian matrix of f{lEv,^)^ calculated in Iev = 1- 

For l,m E {1, 2}, the (Z, m)-th entry of J(l, e), is given by 



J''™(l,e) 



d\ 



EV,l 



^ dlAV,c 



d\ 



EC,e 



diA 



Cm 



(1). (13) 



Consider now a generic {q-^, fc-y) VN of type 7 G ^v having 
at least one codeword socket of edge type e. Using ([8]) and 
(|9]l, it is easy to show that 



lim 

Iav->1 dl 



a 1(7) 

'-' 'ev,; 



(7) 

AV,e 



{lAV,e) 



7,e 



,(7,0l 



,(7,0 



where the last equality is due to 



(1,1) 



we have a, 



(7,1) 



'3'7,i • 



z(i) 



(14) 



0. In fact, for Z = 1 



'^7-i( 10^7,1 ~ 97, i( 10^7,1 ~ ^' analogous way 

0. Next, we develop (fl4] i. assuming 



we can show that a. 



(7,2) 
O.z 



Z = 1, in the two cases Ig = [1,0]"^ and Ig = [0, 1]"^. From 



^ and from the definition of multi-type information function 
in Section Hl-DI we have 



3 z=0 



z=0 
X 



z=0 



S(7) 



|_ 2 '"97,l-2,9-,,2,fe-,-2 

(/c^. - rank(Sg^ i_2. 



)) (15) 



and 



z=0 



[(77,1 97,2 ^7 



^Q-f,l — '^,(l'Y,2 — 'i-,k~f—Z 



E 



(^7 -i'ank(Sg^_i_i,,^ 2 



-l-kj-z)) ■ 



(16) 

It is possible to show that ( fTSl l and (fT6] l are equiv- 
alent to 2 E'Li B^^ill, l)e" = E^=i X^:2(l, l)e" and 
E^=iS^?i(l,2)e" = E^=lX^?^(l,2)e^ respectively. Both 
expressions are obtained through an argument along the same 
line as that used, in the one-edge type case, to prove Lemma 4 
in [7]. Incorporating these expressions into ( fT4] i. recalling that 
\evAIav,£) = J^je^v •^7,ilEV,i(-f^V',e), and recalling 
we finally obtain 



d\EV,l 

OIava 



-(1,^)= E (^)Ex^>,i)^" = ^^^'^(^) 



(17) 



9IeV4 

diAva 



E (^)Ex^>'2KP^"w 



(18) 



Note that in both (fTTl ) and ( fTSl ) the summation is over ^ 



since, for any 7 G J^v \ ^vi, we have X2''^u(lil) ~ 
X2''^(l: 2) = 0. The same proof technique leading to WT\ and 
CI yields §jf^(l,e) ^ P^'He) and §jf^(l,6) = p-'^'{e). 



V2 



We now need to develop q/^"^'" (1) in the right-hand 
side of (l3[ . To this purpose, simply observe that a CN 
of type 6 G may be regarded as a VN whose ks 

local information bits are all punctured (bs = 0). Note that 
this is equivalent to assuming a channel erasure probabil- 
ity e = 1 for the VN. In this case, the right-hand sides 
of m and m become ^Etixitihl) = 



E,e^.2 px'2"i^, 1) and E,.^..^ Er=i 2 
T,.ye.^v2 ^xi^\^,'2), respectively. Thus, we have 

9Iec,i 



dl 



AC,1 



-(1) 



E 



d\ 



EC.l 



(1) 



C,2 



and also ^''^^'^ 



S5i 



er(i,i)=c 



1,1 



1,2 



(19) 



(20) 



(1) = and 



-(1) = C^'i. Hence, 



9Iac,2 ^ ' dlAC,i 

for l,m € {1,2}, the (Z,m)-th entry of J(l,e) is given by 
J''™(1, e) = E;=i P'''(e)C"=''", i.e., J(l, e) = P(e) C □ 



V. Examples 

In this section, the stability of the iterative decoder is 
analyzed for two MET ensembles. 

Example 5.1: Consider the two-edge-type ensemble 
(S' = {1, 2}) whose Tanner graph is depicted in Fig. |2] where 
edges of type 1 ^ S' are depicted in red and edges of type 
2 G in blue. There are N VNs, all of the same type 7. 
Each VN is a length-2 repetition code with = [1, 1], with 
one socket of type 1 e ff and the other of type 2 G ff. 
Thus, we have ~ ^V2 = {7}- There are two CN types 
■!^c — {'5i,(52}, where CNs of type 8\ are (si,/ii) codes, 
depicted in yellow, and CNs of type 82 are [si^h-i) codes, 
depicted in green. All s\ sockets of a type-(5i CN are of type 

1 G (f, while all sockets of a type-(52 CN are of type 2 G (f . 
The number of CNs of types b\ and 62 are N j s\ and N j S2 
respectively, so each edge interleaver is for N edges. 

Assuming that CNs of both types have minimum distance 

2 (^C2 = {(51,^2}), we obtain 



P(6) 



e 
e 



and C 









2 4'^Vs2 



where A'^^' and A^""' are the multiplicities of weight-2 local 
codewords of CNs of types 5i and 82, respectively. From 
Theorem 13.11 the condition for local stability of the erasure- 
free state is 



l(<5: 



e < 



1 



S1S2 



(-51 ) a{S2) 



(21) 



where the right-hand side is an upper bound on the iterative 
decoding threshold called the stability bound. 

From ( l2n i. we see how the multiplicities A^2^'' and ^2''^' 
may jeopardize the decoder stability, and how increasing si 
or S2 is beneficial in terms of stability. We may also observe 
that the erasure-free fixed point for this ensemble is a stable 
attractor if the CNs of at least one type are characterized 




Hi 



Fig. 2. MET ensemble analyzed in Example 15.11 



Fig. 3. MET ensemble analyzed in Example 15.21 



by minimum distance larger than 2, irrespective of the local 
weight spectrum of CNs of the other type (all diagonal entries 
as well as at least one off-diagonal entry of C are zero in 
this case). In practice, for large N this model gives a good 
indication of stability for the ensemble of product codes which 
are obtained by taking N — siS2 and choosing appropriately 
the two edge interleavers Hi and 112 [12]. 

Example 5.2: Consider the two-edge-type ensemble 
{(^ — {1, 2}) whose Tanner graph is depicted in Fig. [3] where 
edges of type 1 G are depicted in red and edges of type 
2 G (o" in blue. There are two VN types = {71 7 72}, 
where the N^-^ VNs of type 71, depicted in cyan, are {q,k) 
codes generated by some generator matrix G^^ , and the N-,,^ = 
Nj-^ q VNs of type 72, depicted in pink, are length-2 repetition 
codes with G^^ — [1,1]. All q sockets of a type-71 VN are 
of type 1 E (a, while both sockets of a type-72 VN are of 
type 2 G Moreover, there are N^-^ q CNs, all of the same 
type 5. Each CN is a (3, 2) SPC code having one socket of 
type 1 G S" and two sockets of type 2 G Hence, we have 

= ^C2 = {S}. 

Assuming that VNs of type 71 have minimum distance 2 

(^V2 = {71 7 72}), we obtain 



P{e) 



2 y^fc 
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and C = 
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1 1 



where B^^u denotes the number of weight- 2 local codewords 
of VNs of type 71 generated by local input words of length 
k through G-y^. Again applying Theorem 13.11 we obtain the 
following condition for stability of the erasure-free fixed point: 



< 1 



(22) 



Again, an increase in the multiplicity of weight-2 local 
codewords of VNs of type 71 has a negative effect on the 
stability of the fixed point / ev ~ 1, as it reduces the range of 
channel erasure probabilities over which such a fixed point is 
locally stable (just note that all coefficients of the polynomial 
on the left-hand side of (l22l) are positive). Moreover, increasing 
q has a positive effect on the stability. We also point out that 
the fixed point Iev = 1 must be locally stable if the minimum 
distance of the type-71 VNs is larger than 2 (in fact, in this 
case we obtain e < 1). Finally, we observe that, upon a proper 
choice of the edge interleaver 112, the Tanner graph depicted in 
Fig-EJcorresponds to the serial concatenation of a (q, k) linear 
block encoder G-., with an accumulator Hence, this class of 



codes may be seen as a generalization of repeat-accumulate 
(RA) codes [13]. An RA code is obtained when type-71 VNs 
are length-g repetition codes. 

VI. Conclusion 

In this paper, the stability condition for iterative BP de- 
coding of MET D-GLDPC codes over the BEC has been 
developed. The obtained inequality is compact, and naturally 
extends to the MET ensemble parametrization the previously 
obtained condition for unstructured irregular (single-edge type) 
D-GLDPC codes. Although this point is not addressed in the 
present work, we mention that for LDPC-like codes, the sta- 
bility condition has a further practical impact on code design 
through its relationship with the average weight distribution 
of the ensemble. 
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